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Abstract. Narrow bandwidth and easy detuning, inefficiency in broadband and 
non-stationary excitations, and difficulties in matching linear harvester’s resonance 
frequency to low-frequency excitations at small scales, have convinced the researchers 
to investigate the nonlinear, in particular the bistable energy harvesters in the recent 
years. However, the bistable harvesters suffer from co-existing low and high energy 
orbits, and sensitivity to initial conditions, and have been recently proven inefficient 
when subjected to many real-world random and non-stationary excitations. Here, we 
propose a novel buy-low-sell-high strategy that can significantly improve the harvester’s 
efficiency at low-frequencies in a much more robust fashion. This strategy could 
be realized by a passive adaptive bistable system. Simulation results confirm high 
efficiency of the adaptive bistable system following a buy-low-sell-high logic when 
subjected to harmonic and random non-stationary walking excitations compared to 
its conventional bistable and linear counterparts. 
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1. Introduction 

Short life span, miniaturization and scalability difficulties, replacement and maintenance 
issues, and relatively very low pace of energy density improvement of conventional 
batteries[1] have convinced many researchers and scientists to consider energy harvesters 
as potential replacement for batteries in many applications. In particular, vibratory 
energy harvesters have captured enormous attention in the last decade due to 
universality and abundant availability of vibratory energy sources. 

Linear harvesters exploiting resonance phenomena suffer from narrow bandwidth 
of efficient harvesting. The narrow resonance bandwidth renders linear harvesters very 
inefficient when subjected to non-stationary excitation where excitation characteristics 
e.g. dominant frequency, change over time, or when the harvester is exposed to 
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broadband random vibration where the excitation power is spread over a wide frequency 
range. To overcome this issue different techniques such as resonance tuning, multi-modal 
energy harvesting, frequency up-conversion, and more recently purposeful inclusion of 
nonlinearity have been suggested [2, 3]. Among these different techniques, deliberate 
introduction of nonlinearity in particular, bistable nonlinearity have been the focus of 
the research in vibration energy harvesting since 2009. ffowever, recent studies have 
revealed that monostable and bistable nonlinear harvesters do not always outperform 
their linear counterparts. 

One of the main issues with the bistable harvester when subjected to harmonic 
excitation is non-uniqueness of the solution and co-existing low-energy and high-energy 
orbits at a given excitation frequency and amplitude [4, 5, 6, 7]. In fact, for a monostable 
nonlinear harvester the probability of converging to the low-energy orbit is higher than 
that of the high-energy orbit [8]. Also, Masana and Daqaq [9] showed that for a given 
excitation level, bistable harvester’s performance is very sensitive to the potential shape 
(shallow versus deep wells). 

Performance of the bistable harvester is further diminished when it is subjected 
to random excitation. Daqaq [10] showed that for an inductive energy harvester with 
negligible inductance, bistability (in general any stiffness nonlinearity) does not provide 
any improvement over the linear one when excited by white noise. Cottone et al. [11] 
and Daqaq [12] showed that when driven by white noise, a necessary condition for 
the bistable harvester to outperform its linear counterpart is to have a small ratio of 
mechanical to electrical time constants. They along with other researchers [13, 14, 15] 
showed that for a given noise intensity, the output power highly depends on the shape 
of the bistable potential. Zhao and Erturk [15] showed that the bistable harvester could 
outperform its linear counterpart only in a narrow region where noise intensity is slightly 
above the threshold of interwell oscillations. 

The bistable harvester becomes even less efficient and less robust when it is excited 
by more realistic and real-world random vibrations (not white noise). Using real 
vibration measurements (human walking motion and bridge vibration) in simulations 
of idealised energy harvesters Green et al. [16] showed that, although the benefits of 
deliberately inducing dynamic nonlinearities into such devices have been shown for the 
case of Gaussian white noise excitations, the same benefits could not be realised for the 
real excitation conditions. 

In this paper, we propose an adaptive bistable harvester that is more robust 
to changes in input excitation parameters and works more efficiently under both 
harmonic and random excitations when compared to its conventional linear and bistable 
counterparts. In the proposed harvester, the potential barrier changes adaptively 
following a buy-low-sell-high strategy [17]. 


Energy harvesting via an adaptive bistable potential 

2. Adaptive bistable harvester 


3 


In this study, we consider both capacitive and inductive harvesters (with single-degree- 


of-freedom in the mechanical domain) with an adaptive bistable potential. Here, the 
adaptive bistable potential refers to a potential where the potential shape, in particular 
the potential barrier height could change according to a logic in an adaptive fashion. 
Adaptive bistability could be realized in different ways that will briefly be discussed 
later in the paper. But first we need to find the logic according to which the bistability 
changes adaptively. The strategy to maximize the harvested energy for a general and 
ideal harvester was derived and discussed in detail by the authors in [17]. However, for 
the sake of readability and its application to the bistable system, its key concepts are 
discussed next. 

2.1. Adaptive bistability logic: Buy-Low-Sell-High 

To find the logic, we consider a model of a single-degree-of-freedom ideal energy harvester 
characterized by the mass m and the displacement x(t) that is subjected to the energy 
harvesting force f(t) and exogenous excitation force F(t). Then, the equation of the 
motion will simply be: 


mx(t ) = F(t) + f{t). 


( 1 ) 


Here we assume that the ideal harvesting force can harvest all the energy that flows 
to the system (there is no accumulation of energy in the system in long-term). Hence, 
maximizing the harvested energy will be equivalent to maximizing the energy flow to 
the system. In other words, we want to maximize: 



( 2 ) 


over admissible trajectories of x(t). It is easy to show that this integral is unbounded if 
x{t) is unconstrained. Indeed, the trajectory defined by a simple relation x{t) = A F{t) 


that can be realized with the harvesting force / = mXF—F results in the harvesting rate 
of A F 2 that can be made arbitrarily large by increasing the mobility constant A. This 
trivial observation illustrates that the question of fundamental limits is only well posed 


for the model that incorporate some technological or physical constraints. This is a 
general observation that applies to most of the known fundamental limits. For example, 
Carnot cycle limits the efficiency of cycles with bounded working fluid temperature, and 


Shannon capacity defines the limits for signals with bounded amplitudes and bandwidth. 


As a common constraint to vibratory energy harvesters, we constrain the harvester 
displacement in a symmetric fashion i.e. |x(t)| < x max , where x max is the displacement 


limit. Rewriting Eq.2 as — f dt F(t)x(t), it could be seen by inspection that the integral 
is maximized by the optimal trajectory: 


x*{t) = —£ max sign F(t) 


( 3 ) 
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This optimal trajectory is indeed realizable by the ideal harvesting force of fit) = 
mx*(t) — Fit). The interpretation of Eq.3 is easy; it says when F(t) is increasing, 
x(t) should be kept at its lowermost limit, and vice versa, when F(t) is decreasing, 
x(t) should be kept at its uppermost limit. Thus, the transitions between displacement 
limits occur when sign of F(t) is changing i.e. at extremums of F(t). In other words, in 
this logic, the harvester mass is kept at its lowest position (—x max ) until the excitation 
force F(t) reaches its maximum when the mass should then be pushed to its highest 
position (x max ) (either by the excitation force, or by the harvesting force if the local 
maximum of the excitation force is still negative or not big enough to push the mass to 
the highest position limit!). Similar dynamics occur in the reverse direction and this 
strategy continues in the same fashion at every extremum of the excitation force F(t). 

If the harvester is incapable of injecting energy to the system (passive-only 
harvester), the harvested mass should traverse between the limits (±x max ) by the 
excitation force F(t) only. In this case, the logic is slightly modified; the harvester 
mass should be kept at its lowest (highest) displacement limit till the largest maximum 
(most-negative minimum) of the excitation force is reached. Only then, the harvester 
mass is pushed from one displacement limit to the other. This logic is very similar to 
the well-known buy-low-sell-high strategy in stock market; hence, we call this logic a 
Buy-Low-Sell-High (BLSH) strategy hereafter. 

Now the question is how to implement this logic. The BLSH strategy could be 
realized by an adaptive bistable potential. In essence, the passive BLSH strategy keeps 
the harvester mass at one end (±x max ) before letting it go to the other end according to 
its logic. A bistable potential with stable points at ±x max and adaptive potential barrier 
could do this. To realize the BLSH logic, the potential barrier should be large enough 
to confine the harvester mass in one well (x max or —x max ). Then, when, according to the 
logic, the harvester mass should travesre to the other end the potential barrier should 
vanish. This logic is schematically shown in Fig.l. 

2.2. Mathematical modeling 

The harvester is modeled as a lumped-parameter mechanical oscillator coupled to a 
simple electrical circuit via an electromechanical coupling mechanism. The formulation 
here is generic and could be applied to both capacitive (e.g. piezoelectric) and inductive 
(e.g. electromagnetic) transduction mechanisms. The nondimensionalized governing 
dynamic equations could be written as [3]: 

dU(x,t) 2 

x + 2(x H-^-- + K 2 y = -x b 

ox 

y + ay = x. (4) 

! It should be noted even though the harvesting force is injecting energy to the system in this case 
during a short period, the net amount of harvested energy will be positive at the end. This is because 
injection of the energy by the harvesting force will pay off when the next excitation force minimum is 
reached. 
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Figure 1 . Passive BLSH strategy realized by an adaptive bistable potential for an 
arbitrary excitation input. The transition from one displacement limit to the other is 
highlighted by the background colour change in the figure. 


In the above equations, x is the oscillator’s displacement relative to base 
displacement ( Xb ). Linear mechanical damping is characterized by the damping ratio 
(, and k denotes the linear electromechanical coupling coefficient, y represents the 
electric quantity that would be voltage or current in capacitive or inductive transduction 
mechanisms, respectively and a is the ratio of the mechanical to electrical time constants. 
The adaptive bistable potential is denoted by U(x, t ) and overdot denotes differentiation 
with respect to nondimensional time. All parameters and variables are nondimensional. 

Two common techniques to realize bistability are buckling phenomenon and 
magnets (to create negative stiffness) in addition to the positive mechanical stiffness. 
When using magnetic field to realize bistability, if permanent magnets are replaced 
by electromagnets [18] (thus having a controllable magnetic field) one can change the 
potential shape; hence, create an adaptive bistability. A passive bistable potential 
admits a quartic form [19], and when made adaptive, we model it as: 

U(x,t) = hi + 5(t)r k )x 2 - j5(t)(l + r k )^, (5) 

Z 4 x s 

where r k < — 1 is strength of the negative stiffness of the magnetic field relative to the 
linear mechanical one. x s denotes the nondimensional stable position of the bistable 
potential, and 5{t) is a Heaviside function which switches between 1 and 0 according to 
the BLSH logic. 5(t) is always equal to unity except when we want the harvester mass 
traverse from one end to the other (according to the BLSH strategy) which then is set 
to zero. 

Figure2(a) depicts an energy harvester with piezoelectric (capacitive) transduction 
mechanism equipped with adaptive bistable potential. The adaptive bistability is 
realized by an electromagnet and a permanent magnet (the proof mass). An On/Off 
controller is used to implement the BLSH logic. The controller senses the excitation and 
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then according to the BLSH strategy sends a signal to the current supplier to supply 
an appropriate current (S(t) = 1) or to shut down the current supply (S(t) = 0). 

Figure2(b) shows how the potential shape changes by the controller signal 5(t), and 
graphically depicts the sequence of the harvester mass trajectory following BLSH logic 
on admissible potential curves§. It should be noted with this type of implementation 
(Eq.5 and Fig. 2(a)) the adaptive bistable system following BLSH logic will not be passive 
for all time. For instance when the harvester mass is moved (T) —> ( 2 ) ((4) —» (3)) a 
positive amount of energy is added to the system because of the way the potential 
shape is changed. However, in the transition right before the one that adds energy, i.e. 
in (2) —* (T) (( 3 ) —* (4)) the same amount of energy is taken out of the system; hence, 
the net energy injected to the system by this type of implementation is zero in half a 
cycle (if not zero for all time) where cycle is referred to transitions from — x max to +x max 
and then back again to — x max . In order to have a passive system for all time, one should 
come up with a bistable mechanism whose potential barrier could be deepened without 
changing the potential energy level of its stable points e.g. like a latching mechanism. 
This is not the case with the current techniques for bistability realization (buckling and 
magnetic field). 

3. Results and discussion 

In this section, simulation results with harmonic and experimental random excitations 
for adaptive bistable harvester is presented and compared with linear and conventional 
bistable harvesters. For a fair comparison, all harvesters are subjected to the same 
displacement limits. To this end, we first optimize the bistable system with respect 
to its potential shape for given excitation input. Then the maximum displacement of 
the optimum bistable harvester is set as the maximum displacement limit for the linear 
and adaptive bistable systems. This approach greatly favors the conventional bistable 
system when it comes to comparison. 

3.1. harmonic excitation 

The potential function considered here for the bistable system is the same as the one 
used for the adaptive bistable harvester with a small change in the parameter notation 
(1 + r k —> —a). The potential used is of the form U(x ) = — \ax 2 + \a ^ where a > 0. 
Fig.3 shows the average power and displacement amplitude of the bistable system when 
subjected to harmonic excitation of the form — Xb = Fosm(uut). This paper intends 
to target mainly the low-frequency excitation where the linear harvesters fail to work 
efficiently; hence, the nondimensional excitation frequency used here is set to tv = 0.05. 
The average power is calculated by \ f*y 2 (t)dt for a long simulation time t. One should 

§ In fact when the magnetic potential is added to the system, the whole bistable potential curve should 
be shifted above the quadratic mechanical potential curve. This does not show up here as we have 
dropped a constant term in Eq.5. However this does not affect the dynamics of the system. 
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(b) 


Figure 2. Energy harvesting with adaptive bistability (a) schematics of a cantilever 
energy harvester with piezoelectric transduction mechanism equipped with adaptive 
bistability (b) change in harvester’s potential function to realize the BLSH logic and 
the sequence of the harvester mass trajectory on admissible potential curves following 
the logic 


note that this expression gives the normalized nondimensional average power. The 
dimensional instantaneous power is equal to {mu^J, 2 )an 2 y 2 where m, co n , and l c are 
the harvester mass, time-scaling frequency, and length scale, respectively. Hence, the 
average power used here is nondimensionalized by mu^l 2 , and further normalized by 

o II 
OLK ||. 

It could be seen from Fig. 3 that the average power increases monotonically with a 
and x s up to a maximum and then drops sharply. This is where the interwell oscillation 
turns into intrawell oscillation (potential barrier linearly increases with a and x 2 ). A 
drastic decrease in the amplitude of oscillation verifies this. It should be noted that for 
values below the optimum value of a (for a given x s ), the system is still in interwell 
motion; however, the power monotonically decreases as a is decreased from its optimum 
value. This could be seen more clearly in Fig. 4. This suggests the robustness issues 
with the conventional bistable system, that is, the harvester works efficiently only when 
the potential barrier is slightly below its critical value when it triggers the interwell 

|| Since we are not optimizing the power with respect to a and n it is fine to normalize the power by 

2 

OLK*. 
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Figure 3. Energy harvesting with conventional bistable system, (a) and (b) show 
surface and contour plots of average harvested power in terms of system parameters a 
and x s . (c) and (d) show surface and contour plots of harvester displacement amplitude 
in terms of system parameters a and x s . The other parameters are set as Fq = 10, 
w = 0.05, C, = 0.01, k = 5, and a = 1000. 



Figure 4. Average harvested power (on the left) and harvester displacement amplitude 
(on the right) of the conventional bistable energy harvester as a function of the potential 
parameter a for three different values of the parameter x s = 2, 3,4. 


oscillation which agrees with Zhao and Etrurk’s claim [15]. 

Next, we compare the performance of the adaptive bistable harvester with that 
of an optimized conventional bistable and linear harvesters when they are subjected to 
harmonic excitation. To this end, we first optimize the parameters of the bistable system 
for given excitation input and displacement limits. The same harmonic excitation used 
in Figs. 3 and 4 is considered here (F 0 = 10 and co = 0.05). According to Figs. 3 and 
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Figure 5. Displacement time histories of linear, conventional bistable, and adaptive 
bistable energy harvesters subjected to harmonic excitation with excitation amplitude 
F o =10, and frequency w=0.05. 


4 the optimal parameters corresponding to maximum displacement of 3.4 are x s = 2 
and a = 12. For a fair comparison the parameters of the adaptive bistable and linear 
harvesters are set such that their maximum displacements do not exceed this value 
[r k = —300 and x s = 2.8 for the adaptive bistable, and natural frequency of \/3 for the 
linear harvester). 

Figures 5 and 6 show time histories of the displacement and electrical-domain state 
(voltage or current for capacitive or inductive transduction mechanisms, respectively) 
for the three adaptive bistable, conventional bistable and linear harvesters. According 
to the figures, although they all have the same maximum displacement, the maximum 
induced voltage (current) in them is quite different with the adaptive bistable having 
the largest and the linear having the smallest induced voltage (current). One could also 
notice the BLSH logic in the adaptive bistable harvester by comparing the moments of 
the transition from one end to the other and the excitation force extremums. It should 
also be noted that the conventional bistable harvester is trying to mimic the BLSH 
strategy in a less effective way. 

Another way to compare the harvesters’ performances is via their phase portraits. 
Fig. 7(a) depicts these phase portraits. As seen in the figure, the transition of the 
oscillator’s mass between the two displacement limits occur at a higher velocity for 
the adaptive bistable harvester than the other two. The force-displacement diagram in 
Fig. 7(b) illustrates it even better as how the adaptive bistable harvester outperforms 
the other two. This diagram shows the force capable of doing positive work versus 
displacement. An ideal harvester i.e. a harvester with BLSH strategy and ideal 
harvesting force, will have a perfect rectangle on this diagram, given the displacement 
limits. This rectangle represents the maximum amount of energy that could be pumped 
into the harvester (which will be consequently harvested by the ideal harvesting force) 
in one cycle. The ideal harvester with the perfect rectangle in the force-displacement 
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Figure 6. Electrical-state (voltage or current depending on transduction mechanism) 
time histories of linear, conventional bistable, and adaptive bistable energy harvesters 
subjected to harmonic excitation with excitation amplitude J 7 o=10, and frequency 
w=0.05. 



Figure 7. Phase portrait (a), and displacement-force diagram (b) of the three 
harvesters when subjected to harmonic excitation with excitation amplitude Fo=10, 
and frequency w=0.05. 


diagram is very analogous to the Carnot cycle with its perfect rectangle in the 
temperature-entropy diagram given the temperature limits of the hot and cold reservoirs. 
In both cases, all the other systems (harvesters and heat engines) fall within this perfect 
rectangle enclosing a smaller area. Time histories of the harvested energy via the three 
harvesters depicted in Fig.8 prove the higher efficiency of the adaptive bistable system 
over the other two. 
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Figure 8. Time history of the harvested energy by the three harvesters when subjected 
to harmonic excitation with excitation amplitude Fo=10, and frequency w=0.05. 


3.2. random excitation: walking motion 

As mentioned earlier, most of the real-world excitations are random and non-stationary 
rather than harmonic, and that the linear and bistable harvesters do not work 
efficiently when subjected to these types of excitations. To examine and compare 
the performance of the three harvesters to random excitations, we subject all the 
harvesters to experimental and relatively low-frequency walking motion. This data 
is experimentally recorded at the hip level while walking [20]. The time history and 
spectral representation of the walking excitation used here are depicted in Fig.9. 

For simulations the experimental data is first non-dimensionalized with scaling 
frequency of 500Hz, and scaling length of 20pm. Again, first the conventional bistable 
potential parameters (a, and x s ) are optimized for maximum harvested energy for a 
displacement constraint of 1.5; then the parameters of the adaptive bistable and linear 
harvesters are set such that they do not exceed this displacement limit. The harvested 
energy is computed the same way as in the case of the harmonic excitation with the 
only difference that it is multiplied by the constant an 2 for the sake of easier numerical 
comparison between different harvesters. 

Fig. 10. (a) illustrates the displacement time history of the harvester with adaptive 
bistability following a BLSH logic. Harvested energy via the harvesters are compared 
in Fig. 10. (b). In addition to the optimal conventional bistable system ( x s = 0.9, and 
a = 1.6), two other bistable systems with detuned a parameter are also simulated. 
According to the figure, BLSH adaptive bistable harvester outperforms the optimal 
conventional bistable and the linear harvesters. It could also be seen that changes in 
the bistable system parameters could significantly diminish the harvester’s efficacy. 







acceleration (m/s 2 ) 
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time (s) frequency (Hz) 


Figure 9. Non-stationary random walking excitation [20]: (a) acceleration time 
history recorded at the hip while walking, and (b) velocity spectrum (Fourier 
transform) of the walking motion 



time 


(b) 



Figure 10. Energy harvesting from walking motion: (a) displacement time history of 
the harvester mass with adaptive bistability subjected to displacement constraint of 
l^maxl < 1-5 (b) energy harvesting time histories of the linear, adaptive bistable, and 
conventional bistable harvesters. Three conventional bistable harvesters with different 
parameters are tested. 




















































































































































Energy harvesting via an adaptive bistable potential 

4. Conclusions 


13 


In this paper, the major drawbacks of the linear and bistable vibration energy harvesters 
were pointed out and a novel adaptive bistable harvester was proposed to overcome 
them. The adaptive bistable harvester follows a buy-low-sell-high strategy. In this 
strategy, the harvester mass is held at the lowest displacement limit (—x max ) before 
the excitation force reaches its maximum which then the harvester mass is pushed 
to the highest displacement limit (x max ) and is held there waiting for the excitation 
force minimum to reach and do the same thing in the reverse fashion. Although this 
strategy guarantees maximum harvested energy in an ideal harvester with no mechanical 
damping, it was shown by simulations that it also works pretty well with more realistic 
set-ups. 

It was also discussed how the adaptive bistable system could be used to enforce 
the BLSH strategy, and how this could be implemented in practice. It was shown that 
a harvester equipped with adaptive bistability following a BLSH logic outperforms its 
linear and conventional bistable counterparts significantly under both harmonic and 
experimental non-stationary random walking excitations. Also the proposed harvester 
does not suffer from the robustness issues similar to those of the linear and conventional 
bistable systems when the system parameters are detuned. Additionally, it was observed 
that at low-frequency excitations the conventional bistable system tries to mimic the 
BLSH strategy which gives an insight to why this harvester is more efficient than its 
linear counterpart at low frequency excitations. 

References 

[1] Anton S R and Sodano H A 2007 Smart materials and Structures 16 R1 

[2] Tang L, Yang Y and Soh C K 2010 Journal of Intelligent Material Systems and Structures 21 

1867- 1897 

[3] Daqaq M F, Masana R, Erturk A and Quinn D D 2014 Applied Mechanics Reviews 66 040801 

[4] Erturk A, Hoffmann J and Inman D 2009 Applied Physics Letters 94 254102 

[5] Mann B and Owens B 2010 Journal of Sound and Vibration 329 1215-1226 

[6] Stanton S C, McGeliee C C and Mann B P 2010 Physica D: Nonlinear Phenomena 239 640-653 

[7] Erturk A and Inman D 2011 Journal of Sound and Vibration 330 2339-2353 

[8] Quinn D D, Triplett A L, Bergman L A and Vakakis A F 2011 Journal of Vibration and Acoustics 

133 011001 

[9] Masana R and Daqaq M F 2011 Journal of Sound and Vibration 330 6036-6052 

[10] Daqaq M F 2011 Journal of Sound and Vibration 330 2554-2564 

[11] Cottone F, Vocca H and Gammaitoni L 2009 Physical Review Letters 102 080601 

[12] Daqaq M F 2012 Nonlinear Dynamics 69 1063-1079 

[13] Litak G, Friswell M and Adhikari S 2010 Applied Physics Letters 96 214103 

[14] Halvorsen E 2013 Physical Review E 87 042129 

[15] Zhao S and Erturk A 2013 Applied Physics Letters 102 103902 

[16] Green P L, Papatheou E and Sims N D 2013 Journal of Intelligent Material Systems and Structures 

24 1494-1505 

[17] Hosseinloo A H and Turitsyn K 2014 arXiv preprint arXiv:lfl0.1f29 

[18] Ouellette S A and Todd M D 2014 Broadband energy harvesting via adaptive control of bistable 



Energy harvesting via an adaptive bistable potential 


14 


potential energy separatrix SPIE Smart Structures and Materials + Nondestructive Evaluation 
and Health Monitoring (International Society for Optics and Photonics) pp 90570Y-90570Y 

[19] Tang L, Yang Y and Soh C K 2012 Journal of Intelligent Material Systems and Structures 

1045389X12443016 

[20] Kluger J M 2014 Nonlinear beam-based vibration energy harvesters and load cells Master’s thesis 



